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^ ' Abstract 
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Fractional analysis is applied to describe classical dynamical systems. Fractional 
derivative can be defined as a fractional power of derivative. The infinitesimal gener- 
ators {H, . } and C = G{q,p)dq + F{q,p)dp, which are used in equations of motion, are 
56 . derivative operators. We consider fractional derivatives on a set of classical observables as 

^H ' fractional powers of derivative operators. As a result, we obtain a fractional generalization 

o ■ 

O 



X 



of the equation of motion. This fractional equation is exactly solved for the simple clas- 
sical systems. The suggested fractional equations generalize a notion of classical systems 
to describe dissipative processes. 
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^ '. 1 Introduction 

^ ■ The analysis of non-integer order goes back to Leibniz, Liouville, Grunwald, Letnikov and 

ly-^ ■ Riemann. There are many books about fractional calculus and fractional differential equations 

[>^ , [D El E] • Derivatives of fractional order, and fractional differential equations have found many 

O I applications in recent studies in physics (see, for example, [HEl ?,[6l ?] and references therein). 

The classical variables, which are also called observables, are defined as functions on the 
phase space. The dynamical description of system is given by an operator. The natural de- 
scription of the motion is in terms of the infinitesimal change of the system. The infinitesimal 
operator of equation of motion is defined by some form of derivation of functions. 

Fractional derivative can be defined as a fractional power of derivative (see, for example, 
[71 [8]). It is known that the infinitesimal generator {H, . } and C = G{q, p)dp + F {q, p)dp, which 
are used in the equation of motion, are derivations on an algebra of classical observables. A 
derivation of an algebra A^ is a linear map C, which satisfies C{AB) = {CA)B + A{CB) for all 
A,B E Ai. In this paper, we consider a fractional derivative as a fractional power of derivative. 
As a result, we obtain a fractional generalization of the equation of motion. It allows us to 
generalize a notion of classical Hamiltonian systems. Note that some fractional generalization 
of gradient systems has been suggested in [9], and a generalization of Hamiltonian systems 
is considered in [10]. The suggested fractional equation is exactly solved for a free particle, 



harmonic oscillator and damped oscillator. A classical system that is presented by fractional 
equation can be considered as a dissipative system. Fractional derivatives can be used as a 
possible approach to describe an interaction between the system and an environment. Note 
that fractional dynamics can be considered with low-level fractionality by some generalization 
of method suggested in [Til [121 [IS] . 

In section 2, the fractional power of derivative and the fractional equation are suggested. In 
section 3, the Cauchy problem for the fractional equation and the properties of time evolution 
described by it are considered. In section 4, the solution for simple examples of the fractional 
equation are derived. 



2 Fractional derivative and fractional equations 

Let us consider the classical systems 

■^qk = Gk{q,p), ■^Pk = Fk{q,p) {k = l,...,n). (1) 

A linear algebra Ai of classical observables is described by functions A = A{q,p) on the phase 
space M^". Let £ be a differential operator on JH given by 

£ = -(G,{q,p)— + F,{q,p)—). (2) 

\ OQk OPk / 

Here and later we mean the sum on the repeated index k from 1 to n. The equation of motion 
for the classical observable has the form 

Equations dTj) are special cases of (jSj). 

If the functions Gk{q,p) and Fk{q,p) satisfy the Helmholtz conditions 

f^-f^^O. (4) 

dpi dpk 

dGi dFk 

TT^ + TT^ = 0' (5) 

oqk opi 

f^-I^^O. ,6) 

oqi oqk 

then the classical system (jH]) is a Hamiltonian system. In this case, Gk and F^ can be represented 

in the form 

^ , . dH OH 



If if = H{q,p) is a continuous differentiable function, then the conditions (jlj), ([5]) and ([6]) are 
satisfied. The equations of motion ([3]) for Hamiltonian system can be written in the form 



j^A = ~{H,A} (7) 



where { , } is the Poisson bracket 



r , „, dA OB dA dB 

The time evolution of the Hamiltonian system is induced by the Hamiltonian H. 

It is interesting to obtain fractional generalizations of equations ([3]) and ([7]). We will consider 
here concept of fractional power for £. If £ is a closed linear operator with an everywhere dense 
domain D{C), having a resolvent R{z,C) = {zLj — C)~^ on the negative half-axis, then there 
exists [m HTJ [18] the operator 



/•oo 

(£)° = ^^^^^ / dzz''~'Ri-z,C)C 

TT Jo 



is defined on D{C) for < a < 1. The operator (£)" is a fractional power of the operator C 
It is known that the linear differential operator ([21) is a closable operator [17]. Note that 



{cncf = (£) 



Q+/3 



for Q!,/3 > 0, q; + /3 < 1. 

As a result, we obtain the equation 

|a = -{crA, (9) 

where t is dimensionless variable. This is the fractional equation of motion. We can define a 
fractional generalization of equation ([7[) by 

|a = -{{H, . }rA, (10) 

Note that {{H, . })" cannot be presented in the form {H', . } with a function H'. Therefore, 
fractional system described by iQ with C = {H, . } are not Hamiltonian systems. The systems 
will be called the fractional Hamiltonian systems (FHS). Usual Hamiltonian systems can be 
considered as a special case of FHS. Note that another fractional generalization of Hamiltonian 
systems has been suggested in [TO] . 



3 Solutions of fractional equations of motion 

3.1 Cauchy problem for fractional equations 

If we consider the Cauchy problem for equation ([3]) in which the initial condition is given at 
the time t = hj Aq, then its solution can be written in the form 

At = $tv4o. 

The operator $i is called the evolution operator. It is not hard to prove that the following 
properties are satisfied: 

<l>t^s = '^t+s, {t,s>0), $0 = /, 

where I A = A for all A & Ai. As a result, the operators $j form a semi-group. Then 
the operator C is called the generating operator, or infinitesimal generator, of the semi-group 
{$i,t>0}. 

Let us consider the Cauchy problem for fractional equation (Q in which the initial condition 
is given by Aq. Then its solution can be presented [171 HB] as 

A{a) = $1")Ao, 

where the operators $J" , t > 0, form a semi-group which will be called the fractional semi- 
group. The operator (£)" is infinitesimal generator of the semi-group {$J" , t > 0} that can be 
presented by 

r(-a) Jo 

This is the Balakrishnan equation [161 113 • 

3.2 Properties of fractional evolution operator 

Let us consider some properties of temporal evolution described by a fractional semi-group 
{$i"\t>0}. 

(1) The operators $J" can be constructed in terms of the operators $< by the Bochner- 
Phillips formula [H [ISl [I?]: 

/■CXD 

*1"^= / dsUt,s)^s, (t>0). (11) 

Jo 

Here fa(t,s) is defined by 

Ut, s) = ^—. dz exp{sz - tz''), (12) 



where a,t > 0, s > 0, and < a < 1. The branch of 2;" is so taken that Re{z") > for 
Re{z) > 0. This branch is a one- valued function in the z-plane cut along the negative real axis. 
The convergence of this integral is obviously in virtue of the convergence factor exp(— tz"). By 
denoting the path of integration in flT2|) to the union of two paths r exp{—i6), and r exp{+i6), 
where r G (0, cxd), and 7r/2 < 6 < it, we can obtain 



1 f^ 

fa(t,s) = — / rfr exp(sr cos^ — tr" cos(q;6'))- 



• sin(sr sin6' — tr"sin(a6') + 6'). (13) 

If we have a solution At of equation ([3]), then formula (TTT]) gives the solution 



Atia)= / ds^(t,s)^, (t>0) (14) 

of fractional equation ([9]). As a result, we can obtain solution of fractional equation by using 
well-known solutions of usual equations. 

(2) In classical mechanics, the most important is the class of real operators. Let * be a 
complex conjugation. If $4 is a real operator on Ai, then 

i^tA)* = <^tiA*) 

for all A G D{^t) G Ai. A classical observable is a real- valued function. If $( is a real 
operator and A is a real-valued function A* = A, then the function At = ^tA is real-valued, 
i.e., {^tA)* = ^tA. An operator, which is a map from a set of observables into itself, should 
be real. All possible dynamics, i.e., temporal evolutions of classical observables, should be 
described by real operators. Therefore the following statement is very important. If $i is a real 
operator, then <l>j is real. The proof will follows from the Bochner-Phillips formula, which 
gives 

POO 

($^)A)*= / dsf:{t,s){^sAy, (t>0). 

Using flT^ . it is easy to see that f*(t, s) = fait, s) is a real- valued function. Then ($iv4)* = ^tA* 
leads to 

i^tA)* = $,(A*) 

for all A G /^(^S"^) C M. 

(3) Let $i be a operator on Ai. An adjoint operator of $i is a operator $4 on Ai*, such 
that 

i^t{A)\B) = (Al^tiB)) 

for all B G -D($i) C Ai and some A E Ai*. The scalar product on At can be defined by 

{A\B)= [ dqdp[A{q,p)]*B{q,p) 



Then an operator $t is called adjoint if 

dqdp[{^tA){q,p)]*B{q,p)= I dqdp[A{q,p)]*{<^tB){q,p). 



Let us give the basic statement regarding the adjoint operator. If $< is an adjoint operator of 
$(, then the operator 

l>J"^= / dsUt,s)^s, (t>0), 



Jo 



is an adjoint operator of $J . We prove this statement by using the Bochner-Phillips formula: 

dqdp ($J"^A)*5 = I ds fait, s) / dqdp ($s^)*5 = 



^0 J 

dsfait,s) f dqdpA*{^sB)= f dqdpA*{<^i''^B). 

The semi-group {<l>t, t > 0} describes a temporal evolution of the distribution function pt{q,p) 
^tPo{Q,p) by the Liouville equation 

Pt{q,p) = -Cpt{q,p), 



dt 
where 



dG,, OF. 



r^r + nfe,rt, o,„„^E(i^-£)- 



If f] < 0, then the system is called dissipative. The semi-group {$J" , t > 0} describes the 
evolution of the density function 

pt{a,q,p) = $l°Vo(g,p) 
by the fractional equation 

This is the fractional Liouville equation. 

(4) It is known that $t is a real operator if $j is real. Analogously, if $j is a real operator, 
then $J° is real. 

(5) Let $j, t > 0, be a positive one-parameter operator, i.e., 

^tA >0 
for A> 0. Using the Bochner-Phillips formula and the property 

fa{t,s)>0 (S>0), 

it is easy to prove that 

$S")a > (A > 0), 

i.e., the operators ^" are also positive. 



4 Examples of fractional equations of motion 

4.1 Fractional free motion of particle 

Let us consider equation ([3]) for free particle. Then 

H = -^/, C = {H, .} = ^d„ 

2m m 

where p is dimensionless variable and m~^ has the action dimension. For A = q, and A = p, 
equation ([3]) gives 

did 

-jTlt = —Pu -irPt = 0. 
dt m dt 

The well-known solutions of these equations are 

t , , 

gt = go H Po, Pt = Po- (15) 

m 

Using these solutions and the Bochner-Phillips formula, we can obtain solutions of the fractional 
equations 

in the form 

poo 

qt{a) = ^i"^qo = / dsU{t,s)qs, Pt{a) = po, 
Jo 

where g^ is given by flT5|) . Then 

gt = go H bait)po, Pt = Po, 

m 

where 

/■oo 
ba{t) = / dsfa{t,s)s. 

Jo 
If a = 1/2, then we have 

'/'^^ 2^Jo V^ 2 

Then 

qt = qo- ^^Poi Pt = Po- (17) 

These equations describe a fractional free motion for a = 1/2. 



4.2 Fractional equation for harmonic oscillator 

Let us consider equation ([3]) for harmonic oscillator. Then L = {H, . }, where 

H = / + q^, (18) 

where t and p are dimensionless variables. For A = q, and A = p, equation ([3]) gives 

-irQt = —Pt, -nPt = -muj'^qt. (19) 

at m at 

The well-known solutions of these equations are 

It = % cos{(jjt) H sin(ci;t), 

Pt = Po cos{ut) — muqo sin^ut). (20) 

Using 02U|) and the Bochner-Phillips formula, we can obtain solutions of the fractional equations 

|g, = -({iJ, . })-g„ |p, = -{{H, . })>, (21) 

where H is defined by f fTSj) . It can be written in the form 

JtP' = ~^ (m^w^gSp - pd^ Pt. 

It is easy to see that these equations with a = 1 give Eqs. (IT9]) . The solutions of fractional 
equations (1^ have the forms 



/■OO 

qt{a) = <l>J°^go = / dsfo,{t,s)qs 
Jo 



.W. 



OO 



Pi(a) = <l>rVo= / dsUt,s)p,. (22) 







Substitution of ([20]) into ([22]) gives 



gt = goCa(t) + — 5a(t), (23) 

Pt = PoCait) - mujqoSa{t), (24) 

where 

/■OO 

Ca('^) = / dsfa{t,s) cos(ci;s), 
Jo 

8 



Sa{t) 



ds fait, s) sm{cos). 



Equations (!23l) and i^M describe solutions of fractional equations fl2T|) for classical harmonic 
oscillator. 

If a = 1/2, then 

^ ,,. t f°° , cosioos) 



Sl/2{t) 



t f^ 



ds 
ds 



-t^/As 



^3/2 

sin(c<;s) 

.s3/2 



-t2/4s 



These functions can be presented through the Macdonald function (see 
such that 



Sec. 2.5.37.1.) 



Cl/2{t) 



/a;t2Ni/4 



V47r 



ujt'^\ 



e+-^/«ir_i/2 ( 2e^-'/\l ^ j + e--^/«ir_i/2 (2e-^"' ^ ^ 



oot' 



Sl/2it) = ^[—^ 



,+ni/8^ 



1/2 



2e--/%/^) - e-/«/^_v2(2e-/y ^ 



where a; > 0, and Ka{z) is the Macdonald function [Ij, which is also called the modified Bessel 
functions of the third kind. 

Note that fractional oscillators are objects of numerous investigations (see, for example, 
[201 ini [221 [231 1211 [25]) because of different apphcations. 



4.3 Fractional equation for damped oscillator 

Let us consider oscillator with linear friction 



d 



d 



dt m dt 

where P < u. The solution of ( !25l) has the form 



-mu'^qt — 2f3p- 



ti 



(25) 



qt = e 



Pt = e 



-pt 



-pt 



go cos 



(Va;2 _ ^2^) + p^sm{^/oj^-(3H) 



muj 



Pocos 



(a/cj^ — j3H) — mujqo sm{^/uj'^ — /3H) 



The fractional equations has the form 



d / V \ '^ 
-T.qt = -{^{mu^q + 2(3p)dp dqj qt, 



d 
Jt 



Pt = -Umuj^q + 2(3p)dp dgj pt. 



(26) 



It is easy to see that these equations with a = 1 give Eqs. (!25l) . Using (1261) and the Bochner- 
Philhps formula, we obtain the solutions 

Pt = PoCa,fiit) - muqoSa,fi(t), (27) 

where 

Jo 

oo 



/■oo 

Jo 



These equations describe solutions of the fractional damped motion of classical harmonic oscil- 
lator. 
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